Enhanced resolution profiling in twins reveals differential methylation signatures of type 2 diabetes with links to its complications.

BACKGROUND: Type 2 diabetes (T2D) susceptibility is influenced by genetic and environmental factors. Previous findings suggest DNA methylation as a potential mechanism in T2D pathogenesis and progression. METHODS: We profiled DNA methylation in 248 blood samples from participants of European ancestry from 7 twin cohorts using a methylation sequencing platform targeting regulatory genomic regions encompassing 2,048,698 CpG sites. FINDINGS: We find and replicate 3 previously unreported T2D differentially methylated CpG positions (T2D-DMPs) at FDR 5% in RGL3, NGB and OTX2, and 20 signals at FDR 25%, of which 14 replicated. Integrating genetic variation and T2D-discordant monozygotic twin analyses, we identify both genetic-based and genetic-independent T2D-DMPs. The signals annotate to genes with established GWAS and EWAS links to T2D and its complications, including blood pressure (RGL3) and eye disease (OTX2). INTERPRETATION: The results help to improve our understanding of T2D disease pathogenesis and progression and may provide biomarkers for its complications. FUNDING: Funding acknowledgements for each cohort can be found in the Supplementary Note.

Body mass index stratified meta-analysis of genome-wide association studies of polycystic ovary syndrome in women of European ancestry.

BACKGROUND: Polycystic ovary syndrome (PCOS) is a complex multifactorial disorder with a substantial genetic component. However, the clinical manifestations of PCOS are heterogeneous with notable differences between lean and obese women, implying a different pathophysiology manifesting in differential body mass index (BMI). We performed a meta-analysis of genome-wide association study (GWAS) data from six well-characterised cohorts, using a case-control study design stratified by BMI, aiming to identify genetic variants associated with lean and overweight/obese PCOS subtypes. RESULTS: The study comprised 254,588 women (5,937 cases and 248,651 controls) from individual studies performed in Australia, Estonia, Finland, the Netherlands and United States of America, and separated according to three BMI stratifications (lean, overweight and obese). Genome-wide association analyses were performed for each stratification within each cohort, with the data for each BMI group meta-analysed using METAL software. Almost half of the total study population (47%, n = 119,584) were of lean BMI (≤ 25 kg/m2). Two genome-wide significant loci were identified for lean PCOS, led by rs12000707 within DENND1A (P = 1.55 × 10-12) and rs2228260 within XBP1 (P = 3.68 × 10-8). One additional locus, LINC02905, was highlighted as significantly associated with lean PCOS through gene-based analyses (P = 1.76 × 10-6). There were no significant loci observed for the overweight or obese sub-strata when analysed separately, however, when these strata were combined, an association signal led by rs569675099 within DENND1A reached genome-wide significance (P = 3.22 × 10-9) and a gene-based association was identified with ERBB4 (P = 1.59 × 10-6). Nineteen of 28 signals identified in previous GWAS, were replicated with consistent allelic effect in the lean stratum. There were less replicated signals in the overweight and obese groups, and only 4 SNPs were replicated in each of the three BMI strata. CONCLUSIONS: Genetic variation at the XBP1, LINC02905 and ERBB4 loci were associated with PCOS within unique BMI strata, while DENND1A demonstrated associations across multiple strata, providing evidence of both distinct and shared genetic features between lean and overweight/obese PCOS-affected women. This study demonstrated that PCOS-affected women with contrasting body weight are not only phenotypically distinct but also show variation in genetic architecture; lean PCOS women typically display elevated gonadotrophin ratios, lower insulin resistance, higher androgen levels, including adrenal androgens, and more favourable lipid profiles. Overall, these findings add to the growing body of evidence supporting a genetic basis for PCOS as well as differences in genetic patterns relevant to PCOS BMI-subtype.

Case Presentation: Functional Assessment of a CASR Variant Identified in a Patient with Hypercalcaemia Confirms Familial Hypocalciuric Hypercalcaemia in the Patient and a Sister Previously Misdiagnosed with Primary Hyperparathyroidism.

Background: Primary hyperparathyroidism (PHPT) and familial hypocalciuric hypercalcaemia (FHH) are common causes of hypercalcaemia. Patients are mostly asymptomatic in the case of FHH and often so in the case of PHPT. In addition, biochemical parameters show considerable overlap, making differential diagnosis difficult. Genetic screening for inactivating variants in the calcium-sensing receptor (CASR) gene that are causative of FHH assists with the diagnosis since such variants are not generally associated with PHPT. However, novel CASR variants must undergo functional assessment before they can be definitively assigned a causative role in FHH. Case Presentations. We describe a 73-year-old female (patient A) who presented with mild parathyroid hormone (PTH)-dependent hypercalcaemia and a history of osteoporosis. Family history revealed that her sister (patient B) had presented a decade earlier with symptoms of PHPT including a history of mild hypercalcaemia and multiple renal calculi, prompting parathyroid surgery. However, a subtotal parathyroidectomy did not resolve her hypercalcaemia long term. On this basis, genetic screening was performed on patient A. This identified a heterozygous variant in the CASR, NM_000388.4:c.T101C: p.Leu34Pro (L34P). Functional analysis showed that the L34P variant was unable to produce mature, dimerized receptor and did not respond to Ca++ ions. Adopting American College of Medical Genetics-based guidelines, the variant was classified as 'Pathogenic (II)'. Patient B was subsequently found to carry the L34P variant heterozygously, confirming a diagnosis of FHH, not PHPT. Conclusion: This study shows the importance of examining patient's family history in providing clues to the diagnosis in isolated cases of hypercalcaemia. In this case, history of a sister's unsuccessful parathyroidectomy prompted genetic screening in a patient who might otherwise have undergone inappropriate parathyroid surgery. Screening detected an inactivating CASR variant, firming up a diagnosis of FHH. These studies reaffirm the requirement for functionally assessing novel CASR variants prior to assigning causality to FHH.

Age-dependent genetic regulation of osteoarthritis: independent effects of immune system genes.

OBJECTIVES: Osteoarthritis (OA) is a joint disease with a heritable component. Genetic loci identified via genome-wide association studies (GWAS) account for an estimated 26.3% of the disease trait variance in humans. Currently, there is no method for predicting the onset or progression of OA. We describe the first use of the Collaborative Cross (CC), a powerful genetic resource, to investigate knee OA in mice, with follow-up targeted multi-omics analysis of homologous regions of the human genome. METHODS: We histologically screened 275 mice for knee OA and conducted quantitative trait locus (QTL) mapping in the complete cohort (> 8 months) and the younger onset sub-cohort (8-12 months). Multi-omic analysis of human genetic datasets was conducted to investigate significant loci. RESULTS: We observed a range of OA phenotypes. QTL mapping identified a genome-wide significant locus on mouse chromosome 19 containing Glis3, the human equivalent of which has been identified as associated with OA in recent GWAS. Mapping the younger onset sub-cohort identified a genome-wide significant locus on chromosome 17. Multi-omic analysis of the homologous region of the human genome (6p21.32) indicated the presence of pleiotropic effects on the expression of the HLA - DPB2 gene and knee OA development risk, potentially mediated through the effects on DNA methylation. CONCLUSIONS: The significant associations at the 6p21.32 locus in human datasets highlight the value of the CC model of spontaneous OA that we have developed and lend support for an immune role in the disease. Our results in mice also add to the accumulating evidence of a role for Glis3 in OA.

Epigenome-wide association study shows differential DNA methylation of MDC1, KLF9 and CUTA in autoimmune thyroid disease.

CONTEXT: Autoimmune thyroid disease (AITD) includes Graves' disease (GD) and Hashimoto's disease (HD), which often run in the same family. AITD etiology is incompletely understood: genetic factors may account for up to 75% of phenotypic variance, whereas epigenetic effects (including DNA methylation (DNAm)) may contribute to the remaining variance (e.g. why some individuals develop GD and others HD). OBJECTIVE: To identify differentially methylated positions (DMPs) and differentially methylated regions (DMRs) comparing GD to HD. METHOD: Whole blood DNAm was measured across the genome using the Infinium MethylationEPIC array in 32 Australian patients with GD and 30 with HD (discovery cohort) and 32 Danish patients with GD and 32 with HD (replication cohort). Linear mixed models were used to test for differences in quantile-normalized beta values of DNAm between GD and HD and data were later meta-analyzed. Comb-p software was used to identify DMRs. RESULTS: We identified epigenome-wide significant differences (p<9E-8) and replicated (p<0.05) 2 DMPs between GD and HD (cg06315208 within MDC1 and cg00049440 within KLF9). We identified and replicated a DMR within CUTA (5 CpGs at 6p21.32). We also identified 64 DMPs and 137 DMRs in the meta-analysis. CONCLUSION: Our study reveals differences in DNAm between GD and HD, which may help explain why some people develop GD and others HD and provide a link to environmental risk factors. Additional research is needed to advance understanding of the role of DNAm in AITD and investigate its prognostic and therapeutic potential.

Leveraging osteoclast genetic regulatory data to identify genes with a role in osteoarthritis.

There has been a growing interest in the role of the subchondral bone and its resident osteoclasts in the progression of osteoarthritis (OA). A recent genome-wide association study (GWAS) identified 100 independent association signals for OA traits. Most of these signals are led by noncoding variants, suggesting that genetic regulatory effects may drive many of the associations. We have generated a unique human osteoclast-like cell-specific expression quantitative trait locus (eQTL) resource for studying the genetics of bone disease. Considering the potential role of osteoclasts in the pathogenesis of OA, we performed an integrative analysis of this dataset with the recently published OA GWAS results. Summary data-based Mendelian randomization (SMR) and colocalization analyses identified 38 genes with a potential role in OA, including some that have been implicated in Mendelian diseases with joint/skeletal abnormalities, such as BICRA, EIF6, CHST3, and FBN2. Several OA GWAS signals demonstrated colocalization with more than one eQTL peak, including at 19q13.32 (hip OA with BCAM, PRKD2, and BICRA eQTL). We also identified a number of eQTL signals colocalizing with more than one OA trait, including FAM53A, GCAT, HMGN1, MGAT4A, RRP7BP, and TRIOBP. An SMR analysis identified 3 loci with evidence of pleiotropic effects on OA-risk and gene expression: LINC01481, CPNE1, and EIF6. Both CPNE1 and EIF6 are located at 20q11.22, a locus harboring 2 other strong OA candidate genes, GDF5 and UQCC1, suggesting the presence of an OA-risk gene cluster. In summary, we have used our osteoclast-specific eQTL dataset to identify genes potentially involved with the pathogenesis of OA.

Genetic determinants of thyroid function in children.

OBJECTIVE: Genome-wide association studies in adults have identified 42 loci associated with thyroid stimulating hormone (TSH) and 21 loci associated with free thyroxine (FT4) concentrations. While biologically plausible, age-dependent effects have not been assessed. We aimed to study the association of previously identified genetic determinants of TSH and FT4 with TSH and FT4 concentrations in newborns and (pre)school children. METHODS: We selected participants from three population-based prospective cohorts with data on genetic variants and thyroid function: Generation R (N = 2169 children, mean age 6 years; N = 2388 neonates, the Netherlands), the Avon Longitudinal Study of Parents and Children (ALSPAC; N = 3382, age 7.5 years, United Kingdom), and the Brisbane Longitudinal Twin Study (BLTS; N = 1680, age 12.1 years, Australia). The association of single nucleotide polymorphisms (SNPs) with TSH and FT4 concentrations was studied with multivariable linear regression models. Weighted polygenic risk scores (PRSs) were defined to combine SNP effects. RESULTS: In childhood, 30/60 SNPs were associated with TSH and 11/31 SNPs with FT4 after multiple testing correction. The effect sizes for AADAT, GLIS3, TM4SF4, and VEGFA were notably larger than in adults. The TSH PRS explained 5.3%-8.4% of the variability in TSH concentrations; the FT4 PRS explained 1.5%-4.2% of the variability in FT4 concentrations. Five TSH SNPs and no FT4 SNPs were associated with thyroid function in neonates. CONCLUSIONS: The effects of many known thyroid function SNPs are already apparent in childhood and some might be notably larger in children as compared to adults. These findings provide new knowledge about genetic regulation of thyroid function in early life.

Bone mineral density loci specific to the skull portray potential pleiotropic effects on craniosynostosis.

Skull bone mineral density (SK-BMD) provides a suitable trait for the discovery of key genes in bone biology, particularly to intramembranous ossification, not captured at other skeletal sites. We perform a genome-wide association meta-analysis (n ~ 43,800) of SK-BMD, identifying 59 loci, collectively explaining 12.5% of the trait variance. Association signals cluster within gene-sets involved in skeletal development and osteoporosis. Among the four novel loci (ZIC1, PRKAR1A, AZIN1/ATP6V1C1, GLRX3), there are factors implicated in intramembranous ossification and as we show, inherent to craniosynostosis processes. Functional follow-up in zebrafish confirms the importance of ZIC1 on cranial suture patterning. Likewise, we observe abnormal cranial bone initiation that culminates in ectopic sutures and reduced BMD in mosaic atp6v1c1 knockouts. Mosaic prkar1a knockouts present asymmetric bone growth and, conversely, elevated BMD. In light of this evidence linking SK-BMD loci to craniofacial abnormalities, our study provides new insight into the pathophysiology, diagnosis and treatment of skeletal diseases.

Digenic Congenital Hypogonadotropic Hypogonadism Due to Heterozygous GNRH1 p.R31C and AMHR2 p.G445_L453del Variants.

A 28-year-old man with congenital hypogonadotropic hypogonadism (CHH) was found to be heterozygous for the GNRH1 p.R31C mutation, reported in the literature as pathogenic and dominant. The same mutation was found in his son at birth, but the testing of the infant at 64 days confirmed the hormonal changes associated with minipuberty. This led to further genetic sequencing of the patient and his son, which found a second variant, AMHR2 p.G445_L453del, in the heterozygous form, reported as pathogenic in the patient but not in his son. This suggests a digenic cause of the patient's CHH. Together, these mutations are postulated to contribute to CHH by the lack of anti-Müllerian hormone (AMH) signalling, leading to the impaired migration of gonadotrophin releasing hormone (GnRH) neurons, the lack of the AMH effect on GnRH secretion, and altered GnRH decapeptide with reduced binding to GnRH receptors. This led us to the conclusion that the observed GNRH1 mutation in the heterozygous state is not certain to be dominant or, at least, exhibits incomplete penetrance and variable expressivity. This report also emphasises the opportunity afforded by the time window of minipuberty in assessing the inherited genetic disorders of hypothalamic function.

Identification of Differentially Expressed Genes and Molecular Pathways Involved in Osteoclastogenesis Using RNA-seq.

Osteoporosis is a disease that is characterised by reduced bone mineral density (BMD) and can be exacerbated by the excessive bone resorption of osteoclasts (OCs). Bioinformatic methods, including functional enrichment and network analysis, can provide information about the underlying molecular mechanisms that participate in the progression of osteoporosis. In this study, we harvested human OC-like cells differentiated in culture and their precursor peripheral blood mononuclear cells (PBMCs) and characterised the transcriptome of the two cell types using RNA-sequencing in order to identify differentially expressed genes. Differential gene expression analysis was performed in RStudio using the edgeR package. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to identify enriched GO terms and signalling pathways, with inter-connected regions characterised using protein-protein interaction analysis. In this study, we identified 3201 differentially expressed genes using a 5% false discovery rate; 1834 genes were upregulated, whereas 1367 genes were downregulated. We confirmed a significant upregulation of several well-established OC genes including CTSK, DCSTAMP, ACP5, MMP9, ITGB3, and ATP6V0D2. The GO analysis suggested that upregulated genes are involved in cell division, cell migration, and cell adhesion, while the KEGG pathway analysis highlighted oxidative phosphorylation, glycolysis and gluconeogenesis, lysosome, and focal adhesion pathways. This study provides new information about changes in gene expression and highlights key biological pathways involved in osteoclastogenesis.

Shared genetics and causal relationships between migraine and thyroid function traits.

BACKGROUND: Epidemiological studies have reported a comorbid relationship between migraine and thyroid dysfunction. METHODS: We investigated the genetic relationship between migraine and thyroid function traits using genome-wide association study (GWAS) data. RESULTS: We found a significant genetic correlation (rg) with migraine for hypothyroidism (rg = 0.0608), secondary hypothyroidism (rg = 0.195), free thyroxine (fT4) (rg = 0.0772), and hyperthyroidism (rg = -0.1046), but not thyroid stimulating hormone (TSH). Pairwise GWAS analysis revealed two shared loci with TSH and 11 shared loci with fT4. Cross-trait GWAS meta-analysis of migraine identified novel genome-wide significant loci: 17 with hypothyroidism, one with hyperthyroidism, five with secondary hypothyroidism, eight with TSH, and 15 with fT4. Of the genes at these loci, six (RERE, TGFB2, APLF, SLC9B1, SGTB, BTBD16; migraine + hypothyroidism), three (GADD45A, PFDN1, RSPH6A; migraine + TSH), and three (SSBP3, BRD3, TEF; migraine + fT4) were significant in our gene-based analysis (pFisher's combined P-value < 2.04 × 10-6). In addition, causal analyses suggested a negative causal relationship between migraine and hyperthyroidism (p = 8.90 × 10-3) and a positive causal relationship between migraine and secondary hypothyroidism (p = 1.30 × 10-3). CONCLUSION: These findings provide strong evidence for genetic correlation and suggest complex causal relationships between migraine and thyroid traits.

Epigenome-Wide Association Study Reveals CpG Sites Associated with Thyroid Function and Regulatory Effects on KLF9.

Background: Thyroid hormones play a key role in differentiation and metabolism and are known regulators of gene expression through both genomic and epigenetic processes including DNA methylation. The aim of this study was to examine associations between thyroid hormones and DNA methylation. Methods: We carried out a fixed-effect meta-analysis of epigenome-wide association study (EWAS) of blood DNA methylation sites from 8 cohorts from the ThyroidOmics Consortium, incorporating up to 7073 participants of both European and African ancestry, implementing a discovery and replication stage. Statistical analyses were conducted using normalized beta CpG values as dependent and log-transformed thyrotropin (TSH), free thyroxine, and free triiodothyronine levels, respectively, as independent variable in a linear model. The replicated findings were correlated with gene expression levels in whole blood and tested for causal influence of TSH and free thyroxine by two-sample Mendelian randomization (MR). Results: Epigenome-wide significant associations (p-value <1.1E-7) of three CpGs for free thyroxine, five for free triiodothyronine, and two for TSH concentrations were discovered and replicated (combined p-values = 1.5E-9 to 4.3E-28). The associations included CpG sites annotated to KLF9 (cg00049440) and DOT1L (cg04173586) that overlap with all three traits, consistent with hypothalamic-pituitary-thyroid axis physiology. Significant associations were also found for CpGs in FKBP5 for free thyroxine, and at CSNK1D/LINCO1970 and LRRC8D for free triiodothyronine. MR analyses supported a causal effect of thyroid status on DNA methylation of KLF9. DNA methylation of cg00049440 in KLF9 was inversely correlated with KLF9 gene expression in blood. The CpG at CSNK1D/LINC01970 overlapped with thyroid hormone receptor alpha binding peaks in liver cells. The total additive heritability of the methylation levels of the six significant CpG sites was between 25% and 57%. Significant methylation QTLs were identified for CpGs at KLF9, FKBP5, LRRC8D, and CSNK1D/LINC01970. Conclusions: We report novel associations between TSH, thyroid hormones, and blood-based DNA methylation. This study advances our understanding of thyroid hormone action particularly related to KLF9 and serves as a proof-of-concept that integrations of EWAS with other -omics data can provide a valuable tool for unraveling thyroid hormone signaling in humans by complementing and feeding classical in vitro and animal studies.

DNA Methylation in Autoimmune Thyroid Disease.

Graves disease and Hashimoto disease form part of the spectrum of autoimmune thyroid disease (AITD), to which genetic and environmental factors are recognized contributors. Epigenetics provides a potential link between environmental influences, gene expression, and thyroid autoimmunity. DNA methylation (DNAm) is the best studied epigenetic process, and global hypomethylation of leukocyte DNA is reported in several autoimmune disorders. This review summarizes the current understanding of DNAm in AITD. Targeted DNAm studies of blood samples from AITD patients have reported differential DNAm in the promoter regions of several genes implicated in AITD, including TNF, IFNG, IL2RA, IL6, ICAM1, and PTPN22. In many cases, however, the findings await replication and are unsupported by functional studies to support causal roles in AITD pathogenesis. Furthermore, thyroid hormones affect DNAm, and in many studies confounding by reverse causation has not been considered. Recent studies have shown that DNAm patterns in candidate genes including ITGA6, PRKAA2, and DAPK1 differ between AITD patients from regions with different iodine status, providing a potential mechanism for associations between iodine and AITD. Research focus in the field is moving from candidate gene studies to an epigenome-wide approach. Genome-wide methylation studies of AITD patients have demonstrated multiple differentially methylated positions, including some in immunoregulatory genes such as NOTCH1, HLA-DRB1, TNF, and ICAM1. Large, epigenome-wide studies are required to elucidate the pathophysiological role of DNAm in AITD, with the potential to provide novel diagnostic and prognostic biomarkers as well as therapeutic targets.

Functional Assessment of Calcium-Sensing Receptor Variants Confirms Familial Hypocalciuric Hypercalcemia.

Context: In the clinic it is important to differentiate primary hyperparathyroidism (PHPT) from the more benign, inherited disorder, familial hypocalciuric hypercalcemia (FHH). Since the conditions may sometimes overlap biochemically, identification of calcium-sensing receptor (CASR) gene variants causative of FHH (but not PHPT) is the most decisive diagnostic aid. When novel variants are identified, bioinformatics and functional assessment are required to establish pathogenicity. Objective: We identified 3 novel CASR transmembrane domain missense variants, Thr699Asn, Arg701Gly, and Thr808Pro, in 3 probands provisionally diagnosed with FHH and examined the variants using bioinformatics and functional analysis. Methods: Bioinformatics assessment utilized wANNOVAR software. For functional characterization, each variant was cloned into a mammalian expression vector; wild-type and variant receptors were transfected into HEK293 cells, and their expression and cellular localization were assessed by Western blotting and confocal immunofluorescence, respectively. Receptor activation in HEK293 cells was determined using an IP-One ELISA assay following stimulation with Ca++ ions. Results: Bioinformatics analysis of the variants was unable to definitively assign pathogenicity. Compared with wild-type receptor, all variants demonstrated impaired expression of mature receptor reaching the cell surface and diminished activation at physiologically relevant Ca++ concentrations. Conclusion: Three CASR missense variants identified in probands provisionally diagnosed with FHH result in receptor inactivation and are therefore likely causative of FHH. Inactivation may be due to inadequate processing/trafficking of mature receptor and/or conformational changes induced by the variants affecting receptor signaling. This study demonstrates the value of functional studies in assessing genetic variants identified in hypercalcemic patients.

Cross-Trait Genetic Analyses Indicate Pleiotropy and Complex Causal Relationships between Headache and Thyroid Function Traits.

Epidemiological studies have reported a comorbid relationship between headache and thyroid traits; however, little is known about the shared genetics and causality that contributes to this association. We investigated the genetic overlap and associations between headache and thyroid function traits using genome-wide association study (GWAS) data. We found a significant genetic correlation (rg) with headache and hypothyroidism (rg = 0.09, p = 2.00 × 10−4), free thyroxine (fT4) (rg = 0.08, p = 5.50 × 10−3), and hyperthyroidism (rg = −0.14, p = 1.80 × 10−3), a near significant genetic correlation with secondary hypothyroidism (rg = 0.20, p = 5.24 × 10−2), but not with thyroid stimulating hormone (TSH). Pairwise-GWAS analysis revealed six, 14, four and five shared (pleiotropic) loci with headache and hypothyroidism, hyperthyroidism, secondary hypothyroidism, and fT4, respectively. Cross-trait GWAS meta-analysis identified novel genome-wide significant loci for headache: five with hypothyroidism, three with secondary hypothyroidism, 12 with TSH, and nine with fT4. Of the genes at these loci, six (FAF1, TMX2-CTNND1, AARSD1, PLCD3, ZNF652, and C20orf203; headache-TSH) and six (HMGB1P45, RPL30P1, ZNF462, TMX2-CTNND1, ITPK1, SECISBP2L; headache-fT4) were significant in our gene-based analysis (pFisher's combined p-value < 2.09 × 10−6). Our causal analysis suggested a positive causal relationship between headache and secondary hypothyroidism (p = 3.64 × 10−4). The results also suggest a positive causal relationship between hypothyroidism and headache (p = 2.45 × 10−3) and a negative causal relationship between hyperthyroidism and headache (p = 1.16 × 10−13). These findings suggest a strong evidence base for a genetic correlation and complex causal relationships between headache and thyroid traits.

Genome-wide analysis of thyroid function in Australian adolescents highlights SERPINA7 and NCOA3.

OBJECTIVE: Genetic factors underpin the narrow intraindividual variability of thyroid function, although precise contributions of environmental vs genetic factors remain uncertain. We sought to clarify the heritability of thyroid function traits and thyroid peroxidase antibody (TPOAb) positivity and identify single nucleotide polymorphisms (SNPs) contributing to the trait variance. METHODS: Heritability of thyroid-stimulating hormone (TSH), free T4 (fT4), free T3 (fT3) and TPOAb in a cohort of 2854 euthyroid, dizygous and monozygous twins (age range 11.9-16.9 years) from the Brisbane Longitudinal Twin Study (BLTS) was assessed using structural equation modelling. A genome-wide analysis was conducted on 2832 of these individuals across 7 522 526 SNPs as well as gene-based association analyses. Replication analysis of the association results was performed in the Raine Study (n = 1115) followed by meta-analysis to maximise power for discovery. RESULTS: Heritability of thyroid function parameters in the BLTS was 70.8% (95% CI: 66.7-74.9%) for TSH, 67.5% (59.8-75.3%) for fT4, 59.7% (54.4-65.0%) for fT3 and 48.8% (40.6-56.9%) for TPOAb. The genome-wide association study (GWAS) in the discovery cohort identified a novel association between rs2026401 upstream of NCOA3 and TPOAb. GWAS meta-analysis found associations between TPOAb and rs445219, also near NCOA3, and fT3 and rs12687280 near SERPINA7. Gene-based association analysis highlighted SERPINA7 for fT3 and NPAS3 for fT4. CONCLUSION: Our findings resolve former contention regarding heritability estimates of thyroid function traits and TPOAb positivity. GWAS and gene-based association analysis identified variants accounting for a component of this heritability.

Familial Dysalbuminemic Hyperthyroxinemia as a Cause for Discordant Thyroid Function Tests.

INTRODUCTION: Discordant thyroid function tests are routinely encountered in clinical practice. Differential diagnoses include acute thyroxine (T4) ingestion, laboratory interference from heterophilic antibodies, thyroid hormone resistance, thyroid-stimulating hormone (TSH)-secreting pituitary adenomas, and T4 protein binding abnormalities. The impact of abnormal binding proteins may be less recognized since widespread use of free T4 (FT4) assays compared to older total T4 assays. CASE REPORT: A 69-year-old female was referred for assessment of discordant thyroid function tests. Biochemistry since July 2015 showed persistently elevated FT4 levels by immunoassay ranging between 25 to 34 pmol/L with normal or slightly decreased TSH ranging between 0.05 to 2.74 mU/L. The patient was clinically euthyroid on 100 mcg daily of levothyroxine for Hashimoto's thyroiditis. FT4 measured using liquid chromatography-tandem mass spectrometry (LC-MS/MS) was 19.5 pmol/L. Exome sequencing (confirmed by Sanger sequencing) detected a guanine to adenine substitution at residue 725 of the ALB gene previously associated with dysalbuminemic hyperthyroxinemia. The patient's daughter had similar thyroid function tests and the same genetic variant. FT4 results from 3 different automated immunoassays showed the Roche Cobas and Siemens Centaur platforms to be most affected by the variant, and Abbott Architect had the best agreement with LC-MS/MS. CONCLUSION: Familial dysalbuminemic hyperthyroxinemia is a potential cause of discordant thyroid function tests. Clinicians suspecting protein-binding abnormalities may further investigate using reference methods such as LC-MS/MS and equilibrium dialysis if available. The increasing accessibility of exome sequencing offers a cost-effective method of diagnosing genetic variants that cause discordant thyroid function tests.

Epigenome-Wide Association Study of Thyroid Function Traits Identifies Novel Associations of fT3 With KLF9 and DOT1L.

CONTEXT: Circulating concentrations of free triiodothyronine (fT3), free thyroxine (fT4), and thyrotropin (TSH) are partly heritable traits. Recent studies have advanced knowledge of their genetic architecture. Epigenetic modifications, such as DNA methylation (DNAm), may be important in pituitary-thyroid axis regulation and action, but data are limited. OBJECTIVE: To identify novel associations between fT3, fT4, and TSH and differentially methylated positions (DMPs) in the genome in subjects from 2 Australian cohorts. METHOD: We performed an epigenome-wide association study (EWAS) of thyroid function parameters and DNAm using participants from: Brisbane Systems Genetics Study (median age 14.2 years, n = 563) and the Raine Study (median age 17.0 years, n = 863). Plasma fT3, fT4, and TSH were measured by immunoassay. DNAm levels in blood were assessed using Illumina HumanMethylation450 BeadChip arrays. Analyses employed generalized linear mixed models to test association between DNAm and thyroid function parameters. Data from the 2 cohorts were meta-analyzed. RESULTS: We identified 2 DMPs with epigenome-wide significant (P < 2.4E-7) associations with TSH and 6 with fT3, including cg00049440 in KLF9 (P = 2.88E-10) and cg04173586 in DOT1L (P = 2.09E-16), both genes known to be induced by fT3. All DMPs had a positive association between DNAm and TSH and a negative association between DNAm and fT3. There were no DMPs significantly associated with fT4. We identified 23 differentially methylated regions associated with fT3, fT4, or TSH. CONCLUSIONS: This study has demonstrated associations between blood-based DNAm and both fT3 and TSH. This may provide insight into mechanisms underlying thyroid hormone action and/or pituitary-thyroid axis function.

Functional Analysis of Calcium-Sensing Receptor Variants Identified in Families Provisionally Diagnosed with Familial Hypocalciuric Hypercalcaemia.

Identification of variants in the calcium-sensing receptor (CASR) gene is an important means of distinguishing between familial hypocalciuric hypercalcaemia (FHH) and primary hyperparathyroidism. However, identification and bioinformatics analysis of genetic variants alone is now considered insufficient as definitive proof; additional functional assessment is required to diagnose FHH with certainty. We identified two novel variants, D433Y and C739Y, and one previously reported variant G509R in the CASR of four kindreds provisionally diagnosed with FHH and aimed to functionally characterise these variants to confirm the diagnosis. Variant receptors were cloned as FLAG-tagged constructs into the mammalian expression vector, pcDNA3.1. Wild type and variant receptor constructs were expressed in HEK293 cells and their expression assessed by Western blot analysis and their functionality analysed using an IP-One assay which measures myo-inositol 1-phosphate accumulation following CaSR activation. Western blot analysis showed that the D433Y receptor had diminished mature glycosylated receptor compared with wild type CaSR whereas the G509R receptor had a complete lack of mature receptor. The C739Y receptor was consistently overexpressed. Functional assessment showed the D433Y receptor to be mildly inactivating at physiological calcium concentrations whereas the G509R receptor was inactive at all calcium concentrations. By contrast, the C739Y variant was activating compared to wild type receptor which is inconsistent with it causing FHH. We conclude that functional assessment of CaSR variants using the IP-One assay was useful in the investigation of suspected FHH probands, confirming the D433Y and G509R variants as likely pathogenic/pathogenic, but dismissing the C739Y variant as causing FHH.